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22 ABSTRACT 

23 Coronaviruses (CoVs) have emerged from animal reservoirs to cause severe and lethal disease in 

24 humans, but there are currently no FDA approved antivirals to treat these infections. One class of 

25 antiviral compounds, nucleoside analogues, mimics naturally occurring nucleosides to inhibit 

26 viral replication. While these compounds have been successful therapeutics for several viral 

27 infections, mutagenic nucleoside analogues, such as ribavirin and 5-fluorouracil, have been 

28 ineffective at inhibiting CoVs. This has been attributed to the proofreading activity of the viral 

29 3’-5’ exoribonuclease (ExoN). |S-D-/V 1 -hydroxycytidinc (NHC, EIDD-1931; Emory Institute for 

30 Drug Development) has recently been reported to inhibit multiple viruses. Here, we demonstrate 

31 that NHC inhibits both murine hepatitis virus (MHV) (EC5o=0.17 uM) and Middle East 

32 respiratory syndrome CoV (MERS-CoV) (ECscf 0.56 pM) with minimal cytotoxicity. NHC 

33 inhibited MHV lacking ExoN-proofreading activity similarly to WT MHV, suggesting an ability 

34 to evade or overcome ExoN activity. NHC inhibited MHV only when added early during 

35 infection, decreased viral specific infectivity, and increased the number and proportion of G:A 

36 and C:U transition mutations present after a single infection. Low-level NHC resistance was 

37 difficult to achieve and was associated with multiple transition mutations across the genome in 

38 both MHV and MERS-CoV. These results point toward a virus-mutagenic mechanism of NHC 

39 inhibition in CoVs and indicate a high genetic barrier to NHC resistance. Together, these data 

40 support further development of NHC for treatment of CoVs and suggest a novel mechanism of 

41 NHC interaction with the CoV replication complex that may shed light on critical aspects of 

42 replication. 
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43 IMPORTANCE 

44 The emergence of corona viruses (CoVs) into human populations from animal reservoirs has 

45 demonstrated their epidemic capability, pandemic potential, and ability to cause severe disease. 

46 However, no antivirals have been approved to treat these infections. Here, we demonstrate the 

47 potent antiviral activity of a broad-spectrum ribonucleoside analogue, P-D-.V '-hydroxycytidi nc 

48 (NHC), against two divergent CoVs. Viral proofreading activity does not markedly impact 

49 sensitivity to NHC inhibition, suggesting a novel interaction between a nucleoside analogue 

50 inhibitor and the CoV replicase. Further, passage in the presence of NHC generates only low- 

51 level resistance, likely due to the accumulation of multiple, potentially deleterious, transition 

52 mutations. Together, these data support a mutagenic mechanism of inhibition by NHC and 

53 further support the development of NHC for treatment of CoV infections. 
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54 INTRODUCTION 

55 The emergence of severe acute respiratory syndrome (SARS) in 2002 and Middle East 

56 respiratory syndrome (MERS) in 2012 has underscored the ability of coronaviruses (CoVs) to 

57 cause lethal disease in humans (1, 2). MERS-CoV continues to infect humans in the Middle East 

58 and four additional human CoVs (HCoVs), HCoV-229E, HCoV-NL63, HCoV-OC43, and 

59 HCoV-HKUl, continue to circulate globally and cause respiratory disease (3-6). The continued 

60 circulation of SARS- and MERS- like CoVs that can replicate efficiently in primary human 

61 airway cells in bat populations further demonstrates the potential for CoVs to emerge and cause 

62 severe disease in the future (7-10). While SARS-CoV and MERS-CoV outbreaks have been 

63 controlled largely through public health measures (11-13), the potential for future outbreaks 

64 highlights the need for safe and effective therapeutics to combat CoV infections. There are 

65 currently no approved therapeutics or vaccines for any human CoV infection. Previous efforts to 

66 treat CoV infections with existing antivirals did not conclusively benefit clinical outcome; thus, 

67 the current standard of care remains mostly supportive (14-16). 

68 Several targets for direct-acting antivirals are being investigated to treat CoV infections 

69 (17-19). Because the viral replication machinery performs an essential role in genome 

70 replication, therapeutics approved to treat multiple different viral infections are aimed at this 

71 target (20). Many approved antivirals are classified as nucleoside analogues, compounds that 

72 mimic natural nucleosides to inhibit viral replication (21). Inhibition by nucleoside analogues 

73 can be accomplished through a variety of mechanisms. Common mechanisms of action include 

74 incorporation of the analogue by the viral polymerase to induce premature termination of strand 
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75 synthesis or loss of essential genetic information through mutagenesis (22-25). A previous study 

76 reported that the nucleoside analogues ribavirin (RBV) and 5-fluorouracil (5-FU) did not 

77 potently inhibit CoVs, and this finding was attributed to the proofreading capabilities of the viral 

78 3’-5’ exoribonuclease (ExoN) (26). Recent reports have demonstrated the inhibition of WT 

79 CoVs by nucleoside analogues such as galidesivir (BCX4430) and remdesivir (GS-5734) (27- 

80 29). While these compounds have shown efficacy against CoVs, administration of multiple 

81 compounds simultaneously may be required to effectively treat CoV infections and control the 

82 emergence of drug resistance, as has been demonstrated for other viral infections (30). 

83 (>-D--hydroxycytidine (NHC) (EIDD-1931, Emory Institute for Drug Development), a 

84 cytidine analogue, has recently been shown to inhibit multiple viruses, including chikungunya 

85 virus, Venezuelan equine encephalitis virus, respiratory syncytial virus, hepatitis C virus, 

86 norovirus, influenza A and B viruses, and Ebola virus (31-36). Previous reports have 

87 demonstrated an increased introduction of transition mutations in viral genomes after treatment 

88 as well as a high genetic barrier to resistance (31, 36). Antiviral activity of NHC has also been 

89 reported against the human a-CoV HCoV-NL63, as well as the (>-CoV SARS-CoV (38-40). 

90 Neither NHC mechanism of action nor NHC resistance have been described for any CoVs to 

91 date. 

92 In this study, we investigated NHC inhibition and resistance in two divergent (S-CoVs, 

93 murine hepatitis virus (MHV) and MERS-CoV. We show that NHC potently inhibits WT MHV 

94 and MERS-CoV with minimal cytotoxicity. We also demonstrate that MHV ExoN proofreading 

95 activity has a limited but measurable effect on sensitivity to NHC. We observe an NHC 

96 inhibition profile consistent with a mutagenic mechanism of action featuring an accumulation of 

97 transition mutations, indicative of a high genetic barrier to resistance. 
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98 RESULTS 

99 NHC inhibits MHV and MERS-CoV replication with minimal cytotoxicity. 

100 NHC (Fig. 1) has potent broad-spectrum antiviral activity against many RNA viral families (31- 

101 36). We first determined if NHC also inhibits CoV replication using a dose-response experiment 

102 with two divergent (3-CoVs: the model CoV, MHV, and the epidemically circulating zoonotic 

103 CoV, MERS-CoV. NHC treatment resulted in a dose-dependent reduction in viral titer for MHV 

104 (Fig. 2A) and MERS-CoV (Fig. 2B). This inhibition resulted in a 50% effective concentration 

105 (EC 50 ) of 0.17 uM for MHV (Fig. 2C) and 0.56 pM for MERS-CoV (Fig. 2D). We detected 

106 negligible changes in DBT-9 cell viability out to 200 pM (Fig. 2E) and CC 50 values above 10 

107 pM in Vero cells (Fig. 2F). The antiviral activity was not due to cytotoxicity, as the selectivity 

108 index was >1000 for MHV and >20 for MERS-CoV. Together, these results confirm potent 

109 inhibition of P-CoVs by NHC. 

110 

111 NHC inhibition profile in CoVs is consistent with mutagenesis. 

112 To better understand the mechanism through which NHC inhibits CoV replication, we performed 

113 a time-of-drug addition assay to determine at what point in the viral replication cycle NHC acts 

114 (37). We added 16 pM (~100x EC 50 concentration) NHC at the indicated times pre- or post- 

115 infection of cells with WT MHV at a MOI of 1 PFU/cell and quantified viral replication after a 

116 single infectious cycle. Compared to the vehicle (DMSO) control, NHC significantly inhibited 

117 MHV replication when added at or before six hours post-infection (Fig. 3A), suggesting that 

118 NHC acts at early stages of the viral replication cycle. We next determined the effect of NHC on 

119 MHV RNA levels and compared to effect on infectious viral titer. RNA levels were reduced by 
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120 approximately 10-fold at the highest tested concentration of NHC in both MHV-infected cell 

121 monolayers (Fig. 3B) and supernatants (Fig. 3C). In contrast, viral titer was reduced up to 5,GOO- 

122 fold at this concentration. We therefore calculated the ratio of infectious virus per viral RNA 

123 genome copy number normalized to the untreated control (specific infectivity) after NHC 

124 treatment and found that the specific infectivity of WT MHV was reduced in a dose-dependent 

125 manner after treatment with increasing concentrations of NHC (Fig. 3D). Together, these data 

126 are consistent with a mutagenic mechanism of NHC anti-CoV activity. 

127 

128 NHC treatment increases transition mutations present across the MHV genome. 

129 To directly test the effect of NHC treatment on mutational burden, we treated WT MHV with 

130 increasing concentrations of NHC and performed full-genome next-generation sequencing 

131 (NGS) on viral populations released after a single round of infection. Our data demonstrate a 

132 dose-dependent increase in mutations present at low frequencies (<5 % of viral population) 

133 across the genome after treatment with increasing concentrations of NHC (Fig. 4A-C). Further 

134 analysis of the types of mutations introduced by NHC revealed an increase in the total number of 

135 transition mutations with increasing NHC concentrations (Fig. 4D-F). The relative proportion of 

136 G:A and C:U transitions among all observed mutations was increased by 13-15% in the presence 

137 of 2 pM NHC and 36-40% in the presence of 4 pM NHC compared to the vehicle control (Fig. 

138 4G, H). Conversely, the relative proportion of A:G and U:C transitions was decreased with 

139 increasing NHC concentrations compared to the vehicle control (Fig. 4G, H). Together, these 

140 results demonstrate that NHC treatment during a single round of WT MHV infection causes 

141 predominantly G:A and C:U transition mutations that are detectable at low frequencies across the 

142 genome. These data further support a mutagenic mechanism of action for NHC inhibition of WT 
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143 MHV. 

144 

145 NHC inhibition is modestly enhanced in the absence of ExoN proofreading. 

146 Mutagenic nucleoside analogues, such as RBV and 5-fluorouracil (5-FU), have been ineffective 

147 at potently inhibiting WT CoVs due to the ExoN proofreading activity (26). A proofreading- 

148 deficient MHV mutant, ExoN(-), displays increased sensitivity to previously tested nucleoside 

149 analogues, indicating that proofreading dampens inhibition by these compounds (26, 38, 39). 

150 Thus, we tested the sensitivity of ExoN(-) MHV to NHC inhibition. Our results indicate that 

151 NHC decreases viral titer of both WT and ExoN(-) MHV in a dose-dependent manner, but that 

152 ExoN(-) MHV demonstrates a statistically significant increase in sensitivity to NHC inhibition 

153 compared to WT MHV (Fig. 5A). However, this difference is reflected in only a modest decrease 

154 in EC 90 concentration by approximately 2-fold for ExoN(-) (0.72 uM) compared to WT MHV 

155 (1.59 pM) (Fig. 5B). The minimal change in sensitivity to NHC observed for ExoN(-) MHV 

156 indicates that NHC potency is only marginally affected by ExoN proofreading activity. 

157 

158 Passage in the presence of NHC yields low-level resistance associated with multiple 

159 transition mutations. 

160 To better understand the development and impact of NHC resistance in CoVs, we passaged two 

161 lineages of WT MHV in thirty times in the presence of increasing concentrations of NHC and 

162 tested the sensitivity of passage 30 (p30) MHV populations to NHC inhibition. We found that the 

163 lineage 1 (MHV p30.1) viral population showed no change in sensitivity to NHC compared to 

164 WT MHV (Fig. 6A). However, lineage 2 (MHV p30.2) showed a decrease in sensitivity to NHC 

165 inhibition in a titer-reduction assay, especially at higher concentrations of compound. We 
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166 observed a modest, approximately 2-fold, increase in EC 90 values for MHV NHC passage 

167 viruses (MHV p30.1 EC 90 = 2.61 pM; MHV p30.2 EC 90 = 2.41 pM; WT MHV EC 90 = 1.53 pM) 

168 (Fig. 6B). This suggests that MHV passage resulted in minimal resistance to NHC. We next 

169 sought to determine if passaging WT MHV in the presence of NHC altered the replication 

170 capacity of these viruses. We found that both lineages showed a delay in replication but 

171 ultimately reached similar peak titers as WT MHV (Fig. 6C). This delay in replication suggests 

172 that MHV p30 viruses are less fit than WT MHV. 

173 To identify mutations associated with these phenotypes after passage, we sequenced 

174 complete genomes of MHV p30.1 and MHV p30.2. Both lineages passaged in the presence NHC 

175 had accumulated over 100 consensus mutations distributed across the genome (Fig. 6D, E; Table 

176 SI). By comparison, a previous study reported that WT MHV accumulated only 23 total 

177 mutations after 250 passages in the absence of drug (39). Further analysis of the p30 MHV 

178 mutational profile demonstrated that slightly more of the total mutations in both lineages were 

179 synonymous changes that did not result in an amino acid change as opposed to nonsynonymous 

180 changes that did alter amino acid sequence (Fig. 6F; Table SI). Additionally, the vast majority of 

181 mutations in both lineages were transition mutations resulting in a purine-to-purine or 

182 pyrimidine-to-pyrimidine change (Fig. 6G). Both lineages contained only two transversion 

183 mutations resulting in a purine-to-pyrimidine or pyrimidine-to-purine change. Though all 

184 possible transition mutation types were detected in both viral lineage populations, the majority in 

185 both passage lineages were G:A transitions (Fig. 6H), which is consistent with the MHV NGS 

186 data (Fig. 4). To determine if the mutational profile at p30 was consistent with an earlier passage, 

187 we analyzed the whole genome of both lineage 1 and 2 at passage 19 (pl9). Both lineages 

188 demonstrated fewer mutations at pl9 than at p30, but the profiles of synonymous vs. 
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189 nonsynonymous changes and the transition mutations were similar (Fig. SI; Table S2). 

190 To determine whether the lack of robust resistance to NHC was broadly applicable across 

191 (3-CoVs, we assessed the capacity of MERS-CoV to evolve resistance to NHC. Like MHV, we 

192 passaged two lineages of MERS-CoV 30 times in the presence of increasing concentrations of 

193 NHC and tested the sensitivity of these lineages to inhibition by NHC. Compared to WT MERS- 

194 CoV passaged in the absence of drug, both MERS-CoV NHC p30.1 and p30.2 exhibited 

195 decreased sensitivity to NHC inhibition (Fig. 7A). This correlated with modestly increased EC 90 

196 values for the passage lineages (WT MERS-CoV EC 9 o= 1.31 pM; MERS-CoV p30.1 EC 9 o = 

197 3.04 pM; MERS-CoV p30.2 EC 9 o= 2.12 pM) (Fig. 7B), corresponding to approximately 2-fold 

198 resistance. Similar to MHV, we observed no substantial shift in dose response curve for MERS- 

199 CoV, indicating minimal acquired resistance. NHC p30 viruses replicated similarly to WT p30 

200 MERS-CoV (Fig. 1C). We sequenced both lineages of MERS-CoV p30 population virus and 

201 detected 27 consensus mutations in MERS-CoV NHC p30.1 (Fig. 7D; Table S3) and 41 

202 consensus mutations in MERS-CoV NHC p30.2 (Fig. 7E; Table S3) that were randomly 

203 distributed across the genome. Both MERS-CoV NHC p30.1 and MERS-CoV NHC p30.2 

204 accumulated nonsynonymous and synonymous mutations in roughly equal proportions (Fig. 7F). 

205 Like in MHV, the mutations detected in MERS-CoV p30 lineages were predominantly transition 

206 mutations (Fig. 7G). Further analysis of these mutations revealed that the predominant type of 

207 transition was lineage-dependent. The majority of transition mutations in MERS-CoV NHC 

208 p30.1 were G:A transitions, as was observed in both p30 MHV lineages, whereas MERS-CoV 

209 NHC p30.2 contained a similar number of each type (Fig. 7H). These results indicate that 

210 MERS-CoV can achieve low-level resistance to NHC and that development of resistance is 

211 associated with the accumulation of multiple transition mutations. Together, our data suggest 
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212 NHC acts as a mutagen and that it poses a high genetic barrier to resistance for fi-CoVs. 

213 DISCUSSION 

214 In this study, we demonstrate that NHC potently inhibits the divergent |J-CoVs MHV and 

215 MERS-CoV. Our data are consistent with a virus mutagenic mechanism of action as evidenced 

216 by a decrease in specific infectivity and an increase in G:A and C:U transition mutations present 

217 at low frequencies across the genome after treatment with NHC. We also demonstrate that robust 

218 resistance to NHC is difficult to achieve in both MHV and MERS-CoV. Both WT MHV and 

219 ExoN(-) MHV are sensitive to NHC inhibition, suggesting that NHC is able to overcome ExoN- 

220 mediated proofreading to inhibit WT CoVs and that it interacts with CoVs differently than other 

221 previously tested nucleoside analogues. 

222 

223 Utility of the broad spectrum antiviral NHC as a pan-CoV therapeutic. 

224 Early work with NHC focused on the mutagenic effects of this compound in multiple bacterial 

225 systems (40-42). More recently, the antiviral properties of this compound have been reported for 

226 multiple RNA viruses, including chikungunya virus, Venezuelan equine encephalitis virus, 

227 respiratory syncytial virus, hepatitis C virus, norovirus, influenza A and B viruses, and Ebola 

228 virus (31-36). NHC has also been shown potently inhibited SARS-CoV and HCoV-NL63 (43, 

229 44), suggesting potential utility in treating CoV infections (17). Based on previous studies, NHC 

230 appears to primarily inhibit viral replication by mutagenesis (31, 34). Serial passaging in the 

231 presence of NHC led to low-level resistance for VEEV, but no detectable resistance for RSV, 

232 IAV, or bovine viral diarrhea virus, indicating a high barrier to resistance (31, 34, 36). Consistent 

233 with these previous studies, we demonstrate that NHC is mutagenic in CoVs and that serial 
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234 passaging yields low-level, approximately 2-fold resistance. Low-level resistance has also been 

235 observed for remdesivir, another nucleoside analogue that potently inhibits CoVs. 

236 Approximately 6-fold resistance to remdesivir is conferred by two mutations in the CoV RdRp 

237 (38). This study further expands the known antiviral spectrum of NHC to include MHV and 

238 MERS-CoV, two genetically divergent (3-CoVs and supports NHC development as a broad- 

239 spectrum CoV antiviral. 

240 

241 NHC inhibition may circumvent ExoN-mediated proofreading. 

242 NHC is the first mutagenic nucleoside analogue demonstrated to potently inhibit proofreading- 

243 intact CoVs. Previous studies have demonstrated that viruses lacking ExoN proofreading 

244 activity, or ExoN(-) viruses, are more sensitive to inhibition by nucleoside analogues, especially 

245 RBV and 5-FU (26, 38, 39, 45). This increased sensitivity has been attributed to an inability of 

246 ExoN(-) to efficiently remove incorrect nucleosides (46). However, we observed a minimal 

247 change in NHC sensitivity between WT MHV and ExoN(-) MHV, especially by EC90. This 

248 suggests that NHC interacts with the CoV replicase differently than other previously tested 

249 nucleoside analogues. One explanation is that NHC may evade removal by the proofreading 

250 ExoN. Studies investigating nucleosides that inhibit DNA viruses have suggested an inability of 

251 the viral exonuclease to efficiently excise some nucleoside analogues (47, 48). Further, a 

252 previous study suggested that the T4 DNA exonuclease activity was incapable of removing NHC 

253 (49). While the SARS-CoV ExoN efficiently removes 3’ terminal mismatches regardless of type 

254 (46, 50), the effect of NHC on this activity has not been investigated. Interestingly, mismatches 

255 readily observed during single nucleotide elongation by the SARS-CoV polymerase in the 

256 absence of drugs correspond to mismatches that would lead to the G:A and C:U transitions 
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257 observed after NHC treatment (46). This suggests that the CoV polymerase could be naturally 

258 more prone to make these types of errors, which are then magnified by NHC. This could lead to 

259 a scenario where ExoN cannot prevent dipping below the error threshold, ultimately resulting in 

260 lethal mutagenesis and similar inhibition of both WT MHV and ExoN(-) MHV (51). 

261 Several nucleosides including the mutagenic RBV have multiple demonstrated 

262 mechanisms other than direct incorporation into the genome (52, 53). Thus, another explanation 

263 for the unique potency of NHC in the presence of an active proofreading ExoN is that it may 

264 inhibit viral replication by additional mechanisms beyond mutagenesis. Indeed, previous reports 

265 have suggested that NHC may also interfere with the RNA secondary structure or virion release 

266 to cause inhibition (31, 36). Further, exogenous C or U in the presence of NHC could rescue 

267 viral replication in in HCV, Chikungunya virus, RSV, and Influenza A virus (32, 34, 36), 

268 indicating that NHC competes with exogenous nucleosides at some stage prior to viral inhibition. 

269 These results raise the possibility that NHC could inhibit a process that results in similar 

270 inhibition of these viruses by a mechanism unrelated to ExoN. Thus, future studies will be 

271 important to investigate the role of proofreading in NHC inhibition of CoVs to shed light on 

272 intricacies of NHC inhibition of the CoV replication complex. 

273 

274 NHC mutagenesis may hinder emergence of robust resistance to NHC. 

275 The decrease in specific infectivity along with the accumulation of transitions across the CoV 

276 genome support a mutagenic mechanism of action for NHC in CoVs. NHC resistance in CoVs 

277 was modest and difficult to achieve, as we obtained approximately 2-fold resistance after 30 

278 passages. Resistance was associated with multiple mutations. Interestingly, MERS-CoV 

279 accumulated less mutations over 30 passages than MHV. While differences in viral mutation 
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280 rates could be the driver of this difference, previous studies have suggested that MHV does not 

281 have a higher mutation rate than MERS-CoV (54-56). The differences in mutation accumulation 

282 between MHV and MERS-CoV may be a product of different passage conditions. While MHV 

283 was passaged with a consistent transfer volume, MERS-CoV passage volumes were adjusted 

284 over time to sustain viral replication under escalating selection for drug resistance. The constant 

285 volume passaging conditions may have more severely bottlenecked MHV populations and fixed 

286 more mutations in the genome than the variable volume passaging conditions applied to MERS- 

287 CoV (57). Alternatively, this difference could also reflect a difference in mutational robustness 

288 of the MHV and MERS-CoV genomes, though this proposition would need to be investigated 

289 further (58, 59). While a portion of the mutations that accumulated over passage likely contribute 

290 to NHC resistance, other mutations, such as those in ns2 or nsp2, which encode proteins 

291 dispensable for viral replication in cell culture, may be merely tolerated because of their limited 

292 effect on viral fitness in the context of our passage conditions (60-62). Few common mutations 

293 arose in both MHV and MERS-CoV passage series, (Supp. Table 1-3), suggesting that multiple 

294 pathways to low-level NHC resistance exist in CoVs. Interestingly, for both MHV and MERS- 

295 CoV, the p30 lineage that demonstrated a greater change in sensitivity to NHC was the lineage 

296 that had fewer overall mutations (Fig. 6, 7). Both MHV passage lineages replicated less well than 

297 WT MHV, suggesting that the accumulation of mutations during passage may negatively impact 

298 viral fitness and the ability of MHV to evolve robust resistance to NHC. Further, the MHV 

299 lineage that did not result in changed sensitivity to NHC by p30 (MHV p30.1) had fewer 

300 mutations present at consensus by pl9 than the other lineage (Fig. SI). Thus, it is possible that 

301 the accumulation of deleterious mutations counteracts potential benefits of resistance mutations 

302 (63). If this is the case, mutations promoting NHC resistance would need to arise early during 
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303 passage to help mitigate the accumulation of excess deleterious mutations. Alternatively, the 

304 inability to evade inhibition by NHC may lead to the accumulation of a greater number of NHC- 

305 associated transitions and ultimately a higher mutational burden that may impact viral fitness 

306 (64, 65). Together, our results support the hypothesis that establishment of resistance to NHC in 

307 CoVs requires a delicate balance of resistance-promoting mutations, viral fitness, and 

308 accumulation of deleterious mutations. Thus, defining the roles of individual NHC resistance- 

309 associated mutations will be an important goal for future studies. Overall, our results support 

310 further development of NHC as a broad-spectrum antiviral for treatment of CoV infections and 

311 contribute new insights into important aspects of CoV replication. 
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312 MATERIALS AND METHODS 

313 Cell culture. Murine astrocytoma delayed brain tumor (DBT) (66) and Vero cells (ATCC CCL- 

314 81) were maintained at 37°C in Dulbecco’s modified Eagle medium (DMEM, Gibco) 

315 supplemented with 10% fetal bovine serum (FBS, Invitrogen), 1% penicillin and streptomycin 

316 (Gibco), and 0.1% amphotericin B (Corning). 

317 Viruses. All work with MHV was performed using the recombinant WT strain MHV-A59 

318 (GenBank accession number AY910861 (67)). MERS-CoV stocks were generated from cDNA 

319 clones (GenBank accession number JX869059(68)). 

320 Compounds and cell viability studies. NHC was synthesized at the Emory Institute for Drug 

321 Development and prepared as a 20 mM stock solution in dimethyl sulfoxide (DMSO). Cell 

322 viability was assessed using CellTiter-Glo (Promega) in 96-well plates according to the 

323 manufacturer’s instructions. DBT and Vero cells were incubated with indicated concentration of 

324 compound at 37°C for 24 hours (DBT) or 48 hours (Vero). Cell viability was determined using a 

325 Veritas Microplate Luminometer (Promega) or GloMax (Promega) with values normalized to 

326 those of vehicle-treated cells. 

327 Nucleoside analogue sensitivity studies and generation of EC S0 curves. Subconfluent 

328 monolayers of DBT cells were infected with MHV at a multiplicity of infection (MOI) of 0.01 

329 PFU per cell for 1 hour at 37°C. The inoculum was removed and replaced with media containing 

330 the indicated compound concentration. Cell supernatants were harvested 24 hours post-infection. 

331 Titers were determined by plaque assay as described previously (69). Subconfluent monolayers 

332 of Vero cells were infected with an MOI of 0.01 PFU/cell of MERS-CoV. After virus adsorption 

333 for 30 minutes at 37°C, the inoculum was removed. Cells were washed with PBS and incubated 

334 with medium containing the indicated concentrations of NHC or DMSO (vehicle control). After 
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335 48 h, supernatant was collected and titers were determined by plaque assay as described 

336 previously (70). EC 50 and EC 90 values and curves were generated using the nonlinear regression 

337 curve fit in GraphPad Prism software (La Jolla, CA). 

338 Time-of-drug addition assay. Subconfluent monolayers of DBT cells were treated with media 

339 containing DMSO or 16 pM NHC (~100x EC 50) ) at the indicated time pre- or post-infection. 

340 Cells were infected with WT MHV at an MOI of 1 PFU/cell for 1 hour at 37°C. Virus inoculum 

341 was removed and fresh medium was replaced. Culture supernatant was harvested 12 hours post- 

342 infection, and viral titer was determined by plaque assay. 

343 Quantification of viral genomic RNA. Subconfluent DBT cells were infected with WT MHV at 

344 an MOI of 0.01 PFU/cell. Inoculum was removed after 1 h incubation at 37°C and medium 

345 containing indicated concentrations of NHC was added. Total RNA from cells and supernatant 

346 RNA was harvested using TRIzol reagent (Invitrogen) after 20 hours. Both total RNA and 

347 supernatant RNA were extracted by phase separation. Total RNA was purified by ethanol 

348 precipitation and supernatant RNA was purified using the PureLink RNA mini kit (Invitrogen) 

349 according to manufacturer’s protocol. Total RNA was reverse transcribed using Superscript III 

350 (Invitrogen) to generate cDNA that was quantified by quantitative polymerase chain reaction 

351 (qPCR) as previously described (26). Data are presented as 2’ AACT ' where AACT denotes the 

352 change in the threshold cycle for the viral target (nsplO) normalized to the control (GAPDH) 

353 before and after drug treatment. Supernatant RNA was quantified using one-step reverse 

354 transcriptase quantitative PCR (qRT-PCR) as previously described (45). Data are presented as 

355 the fold change in genome RNA copies normalized to vehicle control. 

356 Determination of specific infectivity. Subconfluent DBT cells were infected with WT MHV at 

357 a MOI of 0.01 PFU/cell. Inoculum was removed after 1 h incubation at 37°C and medium 
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358 containing indicated concentrations of NHC was added. Supernatant RNA was harvested using 

359 the TRIzol reagent (Invitrogen) after 20 hours, followed by extraction and quantification as 

360 described above. Viral titer was determined by plaque assay. The specific infectivity was 

361 calculated as PFU divided by supernatant genome RNA copy number. This ratio was then 

362 normalized to the vehicle control. 

363 NGS studies. Subconfluent DBT cells were infected with WT MHV at an MOI of 0.01 PFU/cell 

364 and treated with the indicated concentrations of NHC. Supernatant was collected 24 hours post- 

365 infection. Purified viral RNA was submitted to GENEWIZ for library preparation and 

366 sequencing. Briefly, after quality controls, viral RNAs were randomly fragmented using heat. 

367 Libraries were prepared and sequenced on the Illumina HiSeq platform. 

368 GENEWIZ performed base-calling and read demultiplexing. Trimmomatic was used to trim 

369 adapter contaminants and reads shorter than 36 basepairs and filter low quality bases (Q-score 

370 <30) (71). The paired-end fastq reads were then aligned to the MHV genome using Bowtie2 to 

371 generate a SAM file (72). SAMtools was used to process the resultant alignment file and 

372 calculate coverage depth at each nucleotide, generating a sorted and indexed BAM file. LoFreq 

373 was used to call substitution variants, including low-frequency variants, and generate a variant 

374 file (73). The Bash shell and Excel were used to further process and analyze the resultant vcf file. 

375 A frequency of 0.001 was used as a cutoff for variants, consistent with previous reports (74). 

376 Absolute numbers of mutations are reported for each NHC treatment. The percentage of the total 

377 mutations for each specific mutation type was calculated using these numbers. The difference in 

378 percentage for each class of mutation after treatment as compared with vehicle control is referred 

379 to as the relative proportion of these mutations. 

380 MHV population passage in the presence of NHC. WT MHV was passaged in triplicate in 

18 


Downloaded from http://jvi.asm.org/on October 13, 2019 at UNIV OF OTAGO 


cfWology JVI Journal of Virology Accepted MdnUSCfTpf Posted OfTlme 


381 increasing concentrations of NHC, from 1 j_lVT to a maximum of 5 pM. Infection was initiated 

382 for passage 1 at MOI = 0.1 PFU/cell. Viral supernatants were harvested from each viral lineage 

383 and frozen when the cell monolayer demonstrated 80% cytopathic effect (CPE) or after 24 hrs. A 

384 constant volume of 16 pL was used to initiate subsequent passages. All three lineages were 

385 maintained until passage 16 when lineage 3 demonstrated no visible CPE upon multiple attempts 

386 at varying concentrations. Lineage 1 and 2 were maintained until passage 30. After each passage, 

387 total RNA was harvested from infected cell monolayers using the TRIzol reagent. Viral RNA 

388 was extracted from passage 19 and passage 30 samples and reverse transcribed using Superscript 

389 III, followed by generation of 12 PCR amplicons to cover the whole genome. Dideoxy amplicon 

390 sequencing was performed by GENEWIZ and analyzed to identify mutations present at greater 

391 than 50% of total using MacVector. Viral mutation maps depicting the identified mutations were 

392 generated using MacVector. 

393 MERS-CoV population passage in the presence of NHC. Three parallel, independent passage 

394 series of WT MERS-CoV were performed on Vero cells in the presence of gradually increasing 

395 concentrations of NHC up to a maximum concentration of 6.5 pM to select for drug-resistant 

396 mutant viruses. Virus adaptation to NHC-supplemented complete culture medium was assessed 

397 by monitoring progression of characteristic MERS-CoV CPE. Volumes of transferred culture 

398 supernatants were adjusted empirically to balance continuous selective pressure against culture 

399 extinction. Each of triplicate lineages in the MERS-CoV passage experiment was sustained 

400 through passage 30. However, the third lineage was severely impaired in replication and was 

401 excluded from further analysis. Total infected-cell MERS-CoV RNA purified from monolayers 

402 infected with terminal-passage (p30) culture supernatant was used to generate RT-PCR products 

403 for consensus Sanger sequencing of the complete viral genome (Genewiz). Changes in passaged 
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404 virus nucleotide and deduced amino acid sequences were identified via alignment with the WT 

405 parental virus genomic sequence using MacVector. 

406 Virus replication assays. Subconfluent monolayers of DBT (MHV) or Vero (MERS-CoV) cells 

407 were infected with WT or NHC-passaged viral populations at an MOI of 0.01 PFU/cell for 1 

408 hour (MHV) or 30 minutes (MERS-CoV). Inocula were removed and cells were washed with 

409 PBS before addition of pre-warmed media. Supernatants were harvested at indicated times post- 

410 infection, and titers were determined by plaque assay. 

411 Statistics. Statistical tests were performed using GraphPad Prism 7 software (La Jolla, CA) as 

412 described in the respective figure legends. 

413 
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674 Fig. 1. Chemical structure of EIDD-1931, (J - D - A'' 1 - h y cl io x y c y t i cl i n c (NHC). 

675 

676 Fig. 2. NHC inhibits MHV and MERS-CoV with minimal cytotoxicity. (A) Change in MHV 

677 and (B) MERS-CoV titer relative to vehicle control after treatment with increasing 

678 concentrations of NHC. The data represent the results of 6 independent experiments, each with 3 

679 replicates. Error bars represent standard error of the mean (SEM). (C) Change in titer data from 

680 (A) represented as percent of vehicle control. WT MHV EC 5 o= 0.17 pM. (D) Change in titer 

681 data from (B) represented as percent of vehicle control. WT MERS-CoV EC 5 o= 0.56 pM. (E) 

682 DBT-9 cell viability as a percent of DMSO control across NHC concentrations. No cytotoxicity 

683 was detected up to 200 pM. The data represent the results of 2 independent experiments, each 

684 with 2 replicates (MHV). Error bars represent standard error of the mean (SEM). (F) Vero cell 

685 viability as a percent of DMSO control across NHC concentrations. Less than 50% cytotoxicity 

686 was detected up to 10 pM. The data represent the results of 2 independent experiments, each 

687 with 3 replicates. Error bars represent standard error of the mean (SEM). 

688 

689 Fig. 3. NHC inhibition profile of MHV is consistent with mutagenesis. (A) Treatment with 

690 16 pM NHC (~100X EC 50 ) significantly inhibits MHV replication during a single infection when 

691 added before 6 h p.i. (B) Both MHV titer and monolayer RNA copies decrease after treatment 

692 with increasing concentrations of NHC. (C) NHC treatment results in a decrease in supernatant 

693 MHV RNA. (D) Data from (C) represented as the ratio of infectious WT MHV to genomic MHV 

694 RNA present in supernatant, or specific infectivity, normalized to vehicle control. NHC 

695 treatment results in a decrease in specific infectivity of MHV. All data in this figure represent the 

696 results of 2 independent experiments, each with 3 replicates. Error bars represent standard error 

697 of the mean (SEM). Statistical significance compared to DMSO control was determined by one- 
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698 way analysis of variance (ANOVA) with Dunnett’s post-hoc test for multiple comparisons and is 

699 denoted *, P < 0.05; **; P < 0.01; ***, P < 0.001; ****, P < 0.0001. 

700 

701 Fig. 4. NHC treatment drives increase in low-frequency G:A and C:U transition mutations 

702 in WT MHV during a single infection. (A) Distribution and frequency of variants across the 

703 genome detected by NGS after vehicle treatment, (B) 2 pM NHC treatment, or (C) 4 pM NHC 

704 treatment. Logio depth of coverage at each genomic position is depicted by the line; frequency of 

705 individual mutations spread across the genome are represented by dots. (D) Number of mutations 

706 in WT MHV after infection in the presence of (D) vehicle, (E) 2 pM NHC, or (F) 4 pM NHC 

707 presented by type. Transition mutations are shown in grey and transversion mutations are shown 

708 in white. (G) Change in relative proportion of each mutation type after treatment with (G) 2 pM 

709 NHC, or (H) 4 pM NHC compared to vehicle control. The relative proportions of G:A and C:U 

710 transitions increase with increasing concentrations of NHC treatment and are denoted by green 

711 shading. 

712 

713 Fig. 5. Sensitivity of ExoN(-) MHV to inhibition by NHC. (A) Change in viral titer for WT 

714 MHV and ExoN(-) MHV relative to vehicle control after treatment with NHC. ExoN(-) is more 

715 sensitive to NHC than WT. The data represent the results of 3 independent experiments, each 

716 with 3 replicates. Error bars represent SEM. Statistical significance compared to WT MHV was 

717 determined by Wilcoxon test and is denoted **, P < 0.01. (B) Change in viral titer data from (A) 

718 represented as percent of vehicle control. WT ECyo = 1.59 pM, ExoN(-) EC 9 o= 0.72 pM. ExoN(- 

719 ) MHV is approximately 2-fold more sensitive to NHC than WT MHV. 

720 
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721 

722 Fig. 6. Resistance and mutational profile of MHV after 30 passages in the presence of NHC. 

723 (A) Change in viral titer for WT MHV, MHV NHC passage 30 (p30) lineage 1 (MHV p30.1), 

724 and MHV NHC p30 lineage 2 (MHV p30.2) relative to vehicle controls after treatment with 

725 NHC. MHV NHC p30.2 is less sensitive to NHC than WT MHV while MHV p30.1 shows no 

726 change in sensitivity. The data represent the results of 2 independent experiments, each with 3 

727 replicates. Error bars represent SEM. Statistical significance compared to WT MHV was 

728 determined by ratio paired t test and is denoted *, P < 0.05. (B) Change in viral titer data from 

729 (A) represented as percent of vehicle control. WT MHV ECgo= 1.53 pM; MHV p30.1 EC 90 = 

730 2.61 pM, MHV p30.2 EC 9 o= 2.41 pM. (C) Replication kinetics of NHC passage viruses. MHV 

731 p30.1 and p30.2 are delayed in replication compared to WT MHV but ultimately reach similar 

732 peak titers. The data represent the results of 2 independent experiments, each with 3 replicates. 

733 Error bars represent standard deviation (SD). (D) MHV p30.1 accumulated a total of 162 

734 consensus mutations across the genome detectable by Sanger sequencing. Of these mutations, 81 

735 were synonymous, 64 were nonsynonymous, and 17 were noncoding. (E) MHV p30.2 

736 accumulated 102 total mutations across the genome. Of these mutations, 54 were synonymous, 

737 42 were nonsynonymous, and 7 were noncoding. (F) Each lineage accumulated more 

738 synonymous changes than nonsynonymous or noncoding changes over passage. (G) Breakdown 

739 of transition and transversion mutations present in each lineage after passage. MHV p30.1 and 

740 p30.2 mutations were predominantly transitions. (H) Breakdown of the types of transition 

741 mutations present in each lineage across passage. G:A transitions were the most abundant for 

742 both MHV p30.1 and p30.2. 

743 
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744 Fig. 7. Resistance and mutational profile of MERS-CoV after 30 passages in the presence of 

745 NHC. (A) Change in viral titer relative to vehicle controls after treatment with NHC for WT 

746 MERS-CoV passaged 30 times in the absence of drug, MERS-CoV NHC passage 30 lineage 1 

747 (MERS-CoV p30.1), and MERS-CoV NHC passage 30 lineage 2 (MERS-CoV p30.2) relative to 

748 vehicle controls after treatment with NHC. Both MERS-CoV p30.1 and p30.2 are less sensitive 

749 to NHC than WT MERS-CoV. The data represent the results of 2 independent experiments, each 

750 with 3 replicates. Error bars represent SEM. (B) Change in viral titer data from (A) represented 

751 as percent of vehicle control. WT MERS-CoV EC 90 = 1.31 pM, MERS-CoV p30.1 EC 90 = 3.04 

752 pM, MERS-CoV p30.2 EC 9 o= 2.12 pM. (C) Replication kinetics of NHC passage viruses. WT 

753 MERS-CoV, MERS-CoV p30.1, and MERS-CoV p30.2 replicate with similar kinetics and reach 

754 similar peak titers. The data represent the results of 2 independent experiments, each with 3 

755 replicates. Error bars represent SEM.(D) MERS-CoV p30.1 accumulated 27 total mutations 

756 across the genome. Of these mutations, 14 were synonymous and 13 were nonsynonymous. (E) 

757 MERS-CoV p30.2 accumulated 41 total mutations. Of these mutations, 17 were synonymous, 

758 and 24 were nonsynonymous. (F) Both MERS-CoV p30.1 and p30.2 accumulated a similar 

759 numbers of nonsynonymous and synonymous changes during passage. (G) MERS-CoV p30.1 

760 and p30.2 acquired predominantly transitions. (H) The types of transition mutations present in 

761 each lineage across passage. MERS-CoV p30.1 acquired relatively more G:A transitions, 

762 whereas MERS-CoV p30.2 acquired similar numbers of each transition type. 

763 

764 

765 

766 
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